Parametric Forcing of Waves with Non-Monotonic Dispersion 
Relation: Domain Structures in Ferrofluids? 
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Abstract 



Surface waves on ferrofluids exposed to a dc-magnetic field exhibit a non- 
monotonic dispersion relation. The effect of a parametric driving on such 
waves is studied within suitable coupled Ginzburg-Landau equations. Due 



■ to the non-monotonicity the neutral curve for the excitation of standing 

waves can have up to three minima. The stability of the waves with re- 
spect to long-wave perturbations is determined via a phase-diffusion equa- 
tion. It shows that the band of stable wave numbers can split up into two 
or three sub-bands. The resulting competition between the wave numbers 
corresponding to the respective sub-bands leads quite naturally to patterns 
consisting of multiple domains of standing waves which differ in their wave 
number. The coarsening dynamics of such domain structures is addressed. 

February 9, 2008 

05.45.-hb, 05.90.-Fm, 47.20. Tg 



1 



Typeset using REVT^ 



I. INTRODUCTION. 



Spatial patterns have been studied extensively over the past years, the classic examples 
being Rayleigh-Benard convection, Taylor vortex flow, as well as structures arising in di- 
rectional solidification It has been well established that long- wavelength perturbations 
of such steady patterns exhibit diffusive phase dynamics with the phase-diffusion coef- 
ficient changing sign at the Eckhaus instability Thus, after the decay of transients 
one-dimensional steady patterns usually relax to a strictly periodic pattern. 

More recently, it has been pointed out, however, that this need not be the case in general; 
there are situations in which the final state consists of a number of domains with different 
wave numbers This can occur if the phase-diffusion coefficient becomes negative such 

that the band of stable wave numbers is split into two parts. Within each domain the wave 
number is then in one of the two stable sub-bands. Experimentally, such domain structures 



have been observed in Rayleigh-Benard convection in a narrow channel ||T0[- So far it has, 
however, not been firmly established whether the origin of these states is in fact due to a 
splitting of the stable band since the phase-diffusion coefficient has not been measured in 
this regime. From a theoretical point of view this experiment is difficult to analyze since the 
domain structures arise only at very large Rayleigh numbers which require full numerical 
simulations of the three-dimensional Navier-Stokes equations. 

Domain structures can be viewed as arising from the competition between different wave 
numbers. For their investigation it is therefore natural to turn to pattern-forming systems 
which have a neutral curve with two or more minima corresponding to the competing wave 
numbers. Domain structures are then expected to arise close to threshold where they may be 
described analytically within a small-amplitude theory []6|,|9|, p!T|JT^ . In addition, in this regime 
interesting locking |jl3| and coarsening phenomena are possible which are not expected at 
large amplitudes [|n],|12],|14|-|16| . 

So far, not many physical systems have been identified which exhibit a neutral curve 
with two minima. Note that within the framework of Ginzburg-Landau equations it is not 
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sufficient if two different modes with different wave numbers go unstable at the same value 
of the control parameter. Since the stability of the domain structures relies either on phase 
conservation [^,^|12| or on a locking-in of interacting fronts [|1^ (see also below) both minima 
have to correspond to the same mode. 

For traveling waves a neutral curve with a double minimum has recently been identified 
in convection of a conducting fluid in a rotating annulus in the presence of a magnetic field 
[ TTI . This system is, however, not easily accessible in experiments. 

In this paper we show that neutral curves with multiple minima, and the resulting domain 
structures, may be obtained quite naturally by a parametric driving of waves with non- 
monotonic dispersion relation. Then there exists a range of frequencies in which modes with 
different wave numbers resonate simultaneously with the driving. The dispersion relation 
becomes non-monotonic when the group velocity changes sign. This is the case in spiral 



vortex flow between counter-rotating cylinders in certain parameter regimes |T8|. However, 
the parametric forcing of spiral vortex flow has turned out to be non-trivial due to the 
appearance of Stokes layers |]T^ . If the parametric driving can be applied via a bulk forcing 
(e.g. using electric or magnetic fields) its effect is considerably stronger. This suggests 
considering surface waves on ferro-fluids. They exhibit a non-monotonic dispersion relation if 
they are exposed to a sufficiently strong dc- magnetic field In fact, the non-monotonicity 



is a precursor of the Rosensweig instability [^,^. Quite analogously, surface waves on 
conducting and dielectric fluids in the presence of an electric field exhibit a non-monotonic 
dispersion relation leading to an instability for larger fields [p^-p6|. 

The organization of the paper is as follows. To describe parametrically driven waves for 
small amplitudes we introduce in secjl| suitable coupled complex Ginzburg-Landau equa- 
tions. In sec.|T| we study analytically the long-wave stability of standing waves which are 
parametrically excited by a periodic forcing (e.g. by an ac-magnetic or an ac-electric field). 
In the expected parameter regime the band of stable wave numbers separates into 3 sub- 
bands. The numerical simulations presented in sec.|V| con&m the existence and stability of 
structures consisting of domains with different wave number. In addition, we present numeri- 



cal results for the coarsening of arrays of domains within a single, 4* -order Ginzburg-Landau 
equation which can be derived from the coupled Ginzburg-Landau equations under suitable 
conditions. 



II. THE MODEL 

For small amplitudes and small damping parametrically driven waves can be described 
by coupled Ginzburg-Landau equations p7|-|29|l. A crucial ingredient for determining the 
linear part of these equations is the dispersion relation of the waves. For surface waves on 
an inviscid ferro-fluid of infinite depth in the presence of a dc-magnetic field it is given by 
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u^\q) =gq + -q' - n/ ^1/ ^ ^oV- (1) 

Here g is the gravitational acceleration, a the surface tension, p the density of the fluid, 
fi the permeability and Mq the magnetization of the fluid. When Mq is increased beyond 
= y^3^ap( l//io + 1/yu) the dispersion relation becomes non-monotonic and is given by 
a cubic polynomial in the vicinity of the inflection point. To capture this cubic dispersion 
relation third spatial derivatives are retained in the Ginzburg-Landau equations, 

dtA + vd^A = ddlA + fd^^A + aA + hB + cA\A\'^ + {n- c)|5|M, (2) 
dtB - vd^B = d*dlB + f*dlB + a*B + bA + c*B\B\'^ + {n* - c*)\B\^B. (3) 

Physical quantities like the surface height h are described in terms of the complex am- 
plitudes as 

h = 5e^^»^ [A{x, t)e-"'=/2* + B{x, t)e'^^/^') + c.c. + h.o.t., 5 < 1, (4) 

where the amplitudes A and B depend on slow space and time coordinates, t = and 
X = 6x. The parametric driving, which can be achieved with an additional ac-magnetic field, 
enters the equations via the linear coupling terms bA and bB, respectively. Its strength is 
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proportional to the coefficient b which can be chosen real. All other coefficients are 

in general complex. The second control parameter is the detuning between the frequency 
Ue of the external driving and the natural frequency ujq of waves with wave number qq, 
uo = uje/^ + ttj — attr with a being an 0(l)-quantity and a = Qr + iai. The carrier wave 
number go is chosen to be that wave number for which the dispersion relation (|I|) has an 
inflection point with zero slope. This occurs for Mq = Mc- At this point the linear group 
velocity Vr as well as the quadratic dispersion term di vanish. As long as Mq is close to Mc 
these dispersive terms are therefore small and it is consistent to keep also the third-order 
term which gives the cubic dispersion relation. In all of the following we assume that the 
viscosity of the fluid is low. The dissipative terms are then small with the leading order 
term being a^, allowing us to neglect the imaginary part of /. We keep, however, Vi and dr, 
which give the linear and the quadratic dependence of the damping on the wave number, 
although they are also of higher order than a^. In this communication we do not attempt 
to make quantitative predictions for a specific experimental system. Therefore we do not 
calculate the coefficients of (QJ^). Instead we present results that should be typical for the 
parametric driving of waves with non-monotonic dispersion relation. 

The neutral-stability curve, at which the basic state (with flat surface) becomes unstable, 
is given by 

b^ = \Rf with R = a-ivq - dq^ - ifq^. (5) 

The resulting neutral-stability curve is sixth-order in q, and can therefore have up to three 
minima. Depending on the parameters, the absolute minimum of the neutral-stability curve 
can be any one of the three. Of course, the situation with a single minimum can be recovered 
as well. Examples of neutral curves with multiple minima are given by the solid lines in 
figs.|l|(a)-(c). There the effect of changing the driving frequency, i.e. the detuning Oj, is 
demonstrated. The difference between half the external frequency and the frequency of a 
wave with wave number q is indicated by the dotted line. The damping is also taken to 
be wave-number dependent {vi and dr non-zero). As expected the neutral curve exhibits 



local minima at the resonance wave numbers. As the detuning is changed from negative 
to positive values the absolute minimum shifts from the resonance at high wave number to 
that at low wave number. 



III. LINEAR STABILITY AND PHASE DIFFUSION 



The multiple wells in the neutral curve suggest that the band of stable wave numbers also 
splits into separate sub-bands. This is a prerequisite for stable domain structures to exist. 
We therefore determine the linear stability of the nonlinear waves with respect to long-wave 
perturbations. This is done most efficiently by deriving the phase-diffusion equation [|2|fl, 



dT<P{X,T) = V{q)dl<P{X,T). 



(6) 



In this description the amplitudes A and B are proportional to S{q) e^'^l'^ + 0(e) where the 
amplitude 5* of the wave satisfies 



}? = \n\^ S'^ + 2n*RS^ + \R\^ 



(7) 



and X = ex and T = e^t are superslow scales. The phase gives the local wave number via 
q = dx4>- In order to simplify the expression for T>{q) the following notation is introduced 



a * b = ttr-br + ttibi, a • b = + aibr. 
The phase diffusion coefficient 'Dlq) can then be written as D{q)/T{q) with 

D{q) = -{A{3cifq + c*d)\n\'}S'^ 

- {-18c * nV^g^ + 24 (Jm {d*cn)) fq^ 

+ [-12 (Re {v*cn)) f + 8 {d * d*c* n* + 2drdiC • n)] 

+ [6 \nf + 2cin *R) f + 8{Im {cd*nv*))] q 

+2 d * R\nf + 2n* Rc * d — v * v*c * n* — 2vrViC • n ^ 

- {9 [-R *n* + 2RiCi] /2g4 + 12 [ddrRr + {c.di + c*d)Ri + Im (Rnd*)] fq^ 



(8) 



+ [6 {—CiViRr + {ciVr + c • V*) Ri — Re {Rnv*)) f 

+8c * dd * R + Ad * d*n * R* + SdrdiU • R] 
+ [A{v* • Rd*n + v * Rd* •n + c*dv*»R + d* Rv* • c) + QRifn * R]q 
+ 2 {drUrR^ + diUiR^ + d • nRrRij 

-2 (c + n) • RvrV^ + {2ciRi -n*R*) + (2c^i?^ + n * R*) v^] } 
- {9i?2/2g4 ^ * RRifq^ + [qR • + 4{d* Rf] q^ + 4d*RR» v*q +{R» v*f} (9) 

and 

r(g) = -4cr\n\^S^- (AcrH* R + 2Rr\nf) -2Rrn* RS^. (10) 

Here Re{a) and Im{a) denote the real and imaginary part of a, respectively. 

The stability limit (Eckhaus boundary) of the spatially periodic waves is given by 'D^q) = 
and is denoted by the dashed lines in figs.|I](a)-(c). From (|) and (p|,p!OD one can see that 
if the imaginary parts of the various parameters are zero then the neutral-stability and the 
Eckhaus curves are symmetric around q = 0. Even if the neutral curve has multiple minima 
and consequently the band of stable wave numbers is separated into sub-bands, these sub- 
bands merge for large forcing. This can be seen directly from (P|JTI1|) by considering the limit 
S* — > cxD. For Cj = the leading-order term {-4 {3cifq + c*d) \nf] of D{q) is positive 
since Cr < and dr > 0. Thus, T> becomes positive for large S independent of the other 
coefficients, since r > for large S. 

The same consideration shows that for general coefficients, i.e. for q 7^ and / 7^ 0, 
the diffusion coefficient V always becomes negative for large S for some range of q. For the 
parameters chosen in figs.|l|(a)-(c) this occurs for g > 0. Thus, for = —0.7 one obtains two 
stable sub-bands which do not merge for large forcing. With increasing Oj two additional 
resonance arise for smaller q leading to a quite flat neutral curve. Strikingly, the left stable 
sub-band splits up into two bands leading to three separate bands right at onset (flg.|I](b)). 
Upon a further increase to Oj = 0.5 the large-g sub-band merges with the central one for 
small forcing. For large forcing they remain separated. The merging of the central band 



with the left band occurs now only at larger values of the forcing b (fig.|l|(c)). 

IV. NUMERICAL SIMULATIONS 

To show that the complex stability regions shown in figs.|I](a)-(c) indeed lead to stable 
domain structures consisting of domains with large and small wave numbers we have solved 
(@,^ numerically using a finite-difference code with Crank-Nicholson time-stepping. A typi- 
cal domain structure is shown in fig.^. It was obtained for parameter values a = — 0.5 + 0. 5i, 
b = 1.3, c = — 1 — 0.5z, d = 0.1, / = —3, n = —3 + 2i, and v = 1 + 0. 2^. These are the parame- 
ter values that correspond to fig.|l](c) and the location of this domain structure in (g, 6)-space 
is indicated by the dotted line on that figure. The domain structure itself consists of a region 
with q ^ 0.19 in the center of the figure and q ~ —0.72 at the boundaries (thick line). These 
wave numbers correspond to the two Eckhaus-stable wells for the selected parameter values. 
Between the two domains there is a sharp transition. This domain structure is numerically 
stable. The thin line gives Re{exp{iqcx){A + B}) with qc = 1.5 and is intended to give an 
impression of a typical surface deformation in this regime. The broken lines give the real and 
imaginary part of A. The dots in fig.[^(c) give the wave-number distribution of the solution 
shown in fig.^ The higher density of dots near q = —0.7 and q = 0.2 shows that over most 
of the system the local wave number is within one of the two stable regimes. 

From an experimental point of view an important question is how one obtains these 
domain structures. Clearly, an adiabatic increase of the periodic forcing in the frequency 
regime in which the neutral curve has two different minima will not be successful since the 
emerging pattern will always be periodic with the wave number corresponding to the deeper 
minimum. Instead, one has to increase the forcing suddenly from values below threshold to 
a value for which the wave-number band consists of at least two sub-bands. Even then one 
may still predominantly get patterns with the wave number corresponding to the deeper 
minimum. Alternatively, one can change the forcing frequency at fixed supercritical forcing 
amplitude in the regime in which the stable band is split. The change in frequency shifts 
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the band of stable wave numbers and eventually the initially stable wave becomes unstable. 
If its wave number hits the part of the stability boundary which faces the other sub-band 
(as marked by the open square in fig.|l|(c)) it is expected that the instability will not lead to 
a phase slip but to domain structures. 

The best approach is presumably to prepare a periodic pattern in the single-well regime, 
i.e. for M < Mc, with the wave number (i.e. the frequency) chosen such that it falls into 
the unstable region between two sub-bands once the dc-magnetic field is suddenly increased 
to reach M > Mc- The numerical result of such a protocol is shown in fig.|5. It gives the 
space-time diagram for the local wave number and clearly shows how the pattern separates 
into domains with different wave numbers after a jump from Vr > to < which 
renders the initial wave number in the region of instability between the left two sub-bands 
of fig.|I](b) {b = 1.7). In the simulation the initial condition was perturbed with a long- wave 
modulation of the wave number in order to trigger the long-wave Eckhaus instability. Note 
that the initial condition was placed very close to the top of the hump of the Eckhaus curve 
in order to insure that the fastest growing mode has a long wavelength leading to a single 
low-wavenumber domain. Thus, the diffusion coefficient is only weakly negative and the 
evolution is extremely slow. 

A second simulation is shown in fig.§. Here the initial wave number is g = 0.28 and 
6 = 2, i.e. the wave number falls between the two sub-bands which are completely separated 
from each other. The evolution of the wave number is strikingly different from fig.^ Al- 
though the initial perturbation was the same as in fig.|^ large oscillations in the wave number 
arise after a short time. Apparently, in this regime the fastest growing mode has a short 
wavelength. Presumably, this represents the usual shift to shorter wavelengths when the 
Eckhaus boundary is exceeded substantially. In principle it could, however, also indicate an 
additional short-wavelength instability not captured in the phase equation (|^). 

Fig.^ shows that the side-band instability can lead - at least initially - to quite complex 
arrays of domains. As can be seen in that figure, the domains evolve slowly over time in 
a coarsening process in which adjacent domain walls (fronts) annihilate each other thereby 
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reducing the number of domains. Thus, the stabihty of a complex structure hke that shown 
in fig.^ depends crucially on the interaction between the domain walls. Oscillations in the 
wave number - as apparent in fig.0- suggest that the interaction could be oscillatory in space, 
thus allowing for a discrete set of equilibrium distances between the domain walls |13[ . This 
question has been addressed previously within the framework of a single Ginzburg-Landau 
equation with fourth-order spatial derivatives, which models a neutral curve with two equal 



minima 



p!4| ,p!6| , [riip!2| . As discussed there, no locking is possible within the phase equation 
(1). Locking can therefore arise only in regimes in which the phase equation breaks down, 
as is, for instance, the case very close to threshold. 

Close to threshold the coupled Ginzburg-Landau equations (@J^) can be reduced to a 
single, real Ginzburg-Landau equation. To reduce the numerical effort we therefore consider 
now the equation 

dA ^d'^A d^A d^A ^ , , ,,2 , 

^ = D^^ + iDs^ ^ + ^A-\A\A. 11) 

If the neutral curve (|^) is well approximated by a quartic polynomial and if its wells are 
not too deep (|n|) is sufficient to capture the multiple minima In (|TT|) the complex 
amplitude A is proportional to the amplitudes A and B of the left- and right-traveling wave 
components of the standing wave \A\ = \B\ and the cofficients are all real. The slowness 
of the new slow scales X and T is related to the distance from threshold. We expect that 
the results obtained from ([TT|) carry over to the full equations (0,0). In addition, we restrict 



ourselves to the simpler, symmetric case = 0. This case has also been studied in the 
context of two-dimensional zig-zag patterns. There the locking discussed below has been 
identified previously |[l^-[16|. 



In the general, asymmetric case 7^ the minima of the neutral curve have different 
depth. Continuity suggests that the small-amplitude waves near the higher minimum expe- 
rience then a short-wavelength instability which takes them directly to the fastest growing 
mode near the center of the other stability band. For large values of S the stable band 



connected with the higher minimum is known to close ||3^. We therefore expect the results 
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of the symmetric case to carry over to the general case only for intermediate values of S. 

Numerical simulations of eq.(|ll|) yield a surprisingly rich behavior for small domain 
structures which consist only of two domains. This is discussed in detail in ||rT|. Here our 
interest is in the behavior of long systems containing many domains. In particular we are 
interested in the interaction between adjacent domain walls. 

As mentioned above, adjacent domain walls can lock into each other due to their os- 
cillatory interaction in space which is related to the oscillatory behavior of the local wave 
number [|T^ - |l6| , p!l| , |T2| . The solutions of ([TT|) are qualitatively similar to the structure shown 



in fig-lj. The oscillations lead to a discretization of the allowed domain widths, such that 
the domains can be characterized by the number of oscillations contained in them. The 
dependence of the range of existence of a domain on the number of oscillations has been 



studied in [jTT| for domains ranging from 1 to 9 extrema in the wave number. It was found 
that among these domains those with 5 extrema exist over the largest range of D. 

Fig.^(g)-(a) shows the typical evolution of a large array of domains for periodic boundary 
conditions when the coefficient D is increased from negative values. An initial state (fig.^) 
consisting of an arbitrary sequence of domains of different sizes is chosen. Fig.^ shows the 
temporal evolution of the zero crossings of the wave number at fixed D = 0.72; clearly, 
domains that do not have the appropriate width become either wider or narrower in order 
to lock the domain walls at an appropriate distance. The resulting state (fig.^) is stable. In 
analogy to the analysis of spatial chaos in a Ginzburg-Landau equation with real amplitude 



13| , it is expected that in very long states of the type shown in figj^d the sequences of 



domain lengths can be chaotic [|T^,0. 

If D is changed from D = —0.72 to -D = —0.7 states which have only one extremum 
between adjacent walls cease to exist and the adjacent domain walls annihilate each other. 
Such a situation is shown in the sequence fig.^(e)-(c). The single-extremum states disap- 
pear very soon after the control parameter is changed. The remaining states then reorder 
themselves to retain their desired spacing. A similar process is observed in the sequence 
fig.Kc)-(a) when the control parameter is changed to D = —0.58, a value at which the 
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three-extrema states no longer exist. Again, the domain walls annihilate each other and a 
new stable state consisting only of very wide domains is reached. In the numerical simula- 
tions such states with more than 9 extrema were not found to be stable. On an exponentially 
long time scale the state depicted in fig.^ is therefore expected to coarsen to a state with 
only one domain with large and one domain with small wave number. 



V. CONCLUSION 

In this paper we have investigated the stability of parametrically driven waves in sys- 
tems in which the dispersion relation for unforced waves is non-monotonic. We have studied 
coupled Ginzburg-Landau equations which are valid for small amplitudes and in the vicinity 
of the inflection point in the dispersion relation. As expected, the neutral curves for the 
excitation of the waves exhibit multiple minima in this regime, and consequently the band of 
stable wave numbers is split up into three sub-bands. Numerical simulations show that un- 
der these conditions patterns consisting of an array of domains with different wave numbers 
can be stable. Our results suggest that such domain structures should be readily accessible 
experimentally. Particularly suitable appear experiments on surface waves in ferro fluids in 
the presence of a static and a time-periodic magnetic fleld. There the non-monotonic dis- 
persion relation is a precursor of the Rosensweig instability |2l|j2^ . Analogous experiments 
with conducting or dielectric fluids in the presence of the corresponding electric flelds are 
expected to give similar results p3| , ^ , ^ . 

It has been pointed out previously that domain structures are also likely to occur if 
a parametric forcing is applied to waves which are Benjamin- Feir unstable |31[| in the absence 
of forcing. This can be seen from (j^). In the limit of large amplitudes the diffusion coefficient 
becomes negative in the band center if Crdr + Cidi > which is the condition for Benjamin- 



Feir instability of the unforced waves This implies that either the parametrically forced 
waves become unstable at all wave numbers or the stable band has split into subbands. In 
the latter case domain structures should arise. 
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An interesting open question concerns the stability of domain structures in two dimen- 
sions. If the wave numbers in the different domains differ only in their orientation but not in 
their magnitude one obtains zig-zag structures. In isotropic systems they arise generically 
due to the annular shape of the range of stable wave numbers. In axially anisotropic sys- 
tems they appear when the instability of the structureless state is to rolls or waves oblique 
to the preferred direction. The stability of zig-zag structures has been studied in some detail 
rT| , |33[| , in particular in systems with axial anisotropy . 



The neutral stability surfaces for parametrically driven waves in two-dimensional 



anisotropic systems have been studied in some detail in p9[. There interesting situations 
with multiple minima at different orientation as well as different magnitude of the wave num- 
ber have been found suggesting the possibility of a patchwork of two-dimensional domains 
in which the wave numbers differ not only in their orientation but also in their magnitude. 
Can such structures be stable? Some impression can be gained from a simulation presented 
in 



3^ . Fig. 3(b) of [0 shows the result of a two-dimensional simulation of the coupled 
Ginzburg-Landau equations for parametrically driven waves (with normal dispersion) in an 
anisotropic system. The pattern shown is clearly characterized by domains or patches with 
large and small wavelengths. Although the authors attribute the result to a Benjamin- Feir 
instability of the waves, i.e. a situation in which waves of a// wave numbers become unstable, 
comparison with (|^) shows that for the parameters chosen in that simulation the stable band 
is split into two sub-bands with the critical wave number q = lying in the unstable region 
between. Only the stability of the latter wave number had been studied in [|3^. A better 
understanding of such complex two-dimensional domain structures is clearly desirable. 

Finally, an extension to traveling waves is also of interest. So far, domain structures 
of traveling waves have been investigated only within a suitable phase equation The 
possibility of a locking of the domain walls has, however, not been adressed. 
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FIG. 1. The Eckhaus (dashed hne) and neutral-stabihty (sohd hne) curves for eqs.(yj^ for 
ar = —0.5, d = 0.1, V = 1 + 0.2i, c = —1 — 0.5i, n = —3 + 2i and various values of the detuning 
as indicated. The dotted hne in (c) indicates the wavenumber distribution of the structure shown 
in figJI 
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FIG. 2. A typical domain structure arising from eqs. (IH) for ar = -0.5 + 0.5i, d = 0.1, 
V = 1 + 0.2i, c = -1 - 0.5i, n = -3 + 2i. The thin solid line gives i?e(e*''=''(A + B)) with = 1.5. 
The short and long dashed lines indicate the real and imaginary part of A, respectively. 
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FIG. 3. Temporal evolution of the local wave number. The initial wave number is q m —0.34. 
The forcing is 6 = 1.7. The remaining parameters are a = —0.5 + O.li, d = 0.1, v = 1 + 0.2i, 
c = -1 - 0.5i, n = -3 + 2i, f = -3. 
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FIG. 4. Temporal evolution of the local wave number. The initial wave number is q = 0.28 and 
6 = 2. The remaining parameters are as in fig.^ The fastest growing mode has a short wavelength 
and leads to a large number of small domains which then undergo a coarsening process. 
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FIG. 5. Positions of domain walls as they evolve, and the domain structures observed at various 
times. Parts (a), (c), (e), and (g) show the wave number (thick line) and Re{Ae^'^°^) (thin line) 
while (b),(d), and (f) plot the positions of zeros of the wave number. For all figures, S = 1.0 and 
for (a-c) D = -0.72, (d,e) D = -0.7 and (f,g) D = -0.58. The times indicated in figures (b), (d), 
and (f) are time/10, 000. 
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